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Introduction. In this work we investigate the cascade fragmentation 
process of material sample. This process is the particular case of the cascade 
kinetic processes which arise in different fields of theoretical physics (see, for 
example, [I]). In spite of the fact that cascade type of fragmentation is the 
simplest one and the fragmentation which has no the cascade character is 
more widespread in nature, such processes are important from the physical 
point of view and have been not completely studied. In particular, it is not 
studied the cascade fragmentation process that is not completely self-similar 
in contrast to the analyzed one in the classical work [2]. In the present 
communication, we investigate the fragmentation being not self-similar that 
J±« it is analogous to the process introduced in the paper [3j. 

General equation. We describe the material fragmentation process by 
the particle density n(r,t). The value n(r,t)dr is the average fragment num- 
ber which have their sizes in the interval (r, r + dr) at the time moment t. 
The fragmentation process is determined completely by the conditional prob- 
ability density P(p,r,i) of the splinter formation of given fragment. More 
strictly, the probability of the formation of splinter that has the size in the 
interval (p, p + dp) during the time interval (t,t + dt) is equal to P(p, r, t)dpdt 
when it is occurred the decay of the fragment with the size r. We introduce 
also the intensity p(r,t) of the fragment decay. It is the average number of 
decays of fragments with the size r in the time interval (t, t + dt). 

The evolution of the density n(r, t) is described in terms of values pointed 



out and on the base of the kinetic equation 

oo 

n(r,t) = Jp(r, P ,tH P ,t)d P -Kr,tHr,t). (1) 

r 

In the article [3] , it has been studied the fragmentation process with the scale- 
invariant subdivision having temporally independent conditional probability 
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density of transitions 

P(p,r,t) = -P(p/r). 
r 

It has been done following the idea contained in the pioneer work [2] . At the 
present communication, it is considered the fragmentation process that has 
no the scale-invariant subdivision when the transition probability density has 
the form 

P(p, r,t) = P(p/r) (2) 

as it has been done in the paper [3j. So, the chosen subdivision scale exists 
since, in this case, the density P(p, r, t) has the physical dimensionality of the 
inverse length. Following to works pf], [5], we take into account the volume 
conservation in the fragment system 

DC 

V = J r 3 n(r, t)dr = const . (3) 
o 

This property gives the next relation for the intensity 

r 1 

p(r,t) = J P(p/r) ^-^ dp = r J P(x)x 3 dx = pr, p = const. (4) 



Thus, we consider the next kinetic equation for the particle density 
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n(r,i) = / P(r / p)n(p,t)dp — prn(r,t) . (5) 



The limit distribution. The formula (2) permits to use the Mellin 
transformation 

00 

.s-l„ 



M(s,t) = / r s ' l n{r,t)dr 



o 



for study some solutions of the equation (5). In this case, M(l,t) = N(t) 
is the average fragment number at the moment t; M(2,t) = N(t)X(t) where 
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X(t) is the average fragment size; M(3, t) = S(t) is the total fragment surface 
M(4, t) = V = const is the total fragment volume. 
The application of the Mellin transformation to both parts of the equation 
(5) leads to the equation 

M(s,t) = M(s + l,t)(p(s)-n), (6) 

1 

where p(s) = J x s ~ 1 P(x)dx and \i = p(4) according to (4). In particular, 
o 

M(3, t) = M(3, 0) + V(p(3) - p(4))t (7) 

since M(4, t) = V = const. The asymptotic behavior of the distribution 
f(r, t) density at t — > oo is of interest for us. In this case, it is sufficient to 
use the second term M(3, t) ~ V(p(3) —p(4))t only in the last formula. From 
equations (6) and (7), we find 

M(2,*) - iy(p(3) -p(4))(p(2) -p(4))t 2 = N(t)\(t), (8) 

M(l,t) - iy(p(3) -p(4))(p(2) -p(4))(p(l) -p(4))* 3 = . (9) 
Therefore, it follows 

Af(2,Q 3 1 

AW ~M(M) = ?aF^(4)*' (10) 

The analysis of the model P(x) ~ x a . In general case, the study of 
solutions of the differential equation (6) is sufficiently difficult. However, for 
the probability distribution density f(r,t) = n(r,t)/N(t) on fragment sizes 
when the variable r is replaced by X(t)r, the analysis of this equation becomes 
much more simple. Let us introduce the designation by F(s,t) of the Mellin 
transformation 

oo 



F(s,t) = \(t) J r^f^Xity^di 
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of the distribution X(t)f(X(t)r,t) density. Changing the variables r =4> X(t)r 
in the equation (5) and computing the Mellin transformation pointed out, we 
obtain 

d /JV\ s 

F(s, t) = -F(s, t)-\n l-\ - -F(s, t)X + X(p(s) - ^)F(s + l,t). (11) 

The distribution X(t)f(X(i)r,t) density and its Mellin transformation F(s,t) 
are tended to /oo(r) and F(s) at t — > oo, correspondingly. These limit func- 
tions are nontrivial 0;oo). In this case, F(s,t) — > 0, £ — > oo. So, the 
equation (11) take the form 

= -F(s)± In (£) + A(p( S ) - M )F( S + 1) - A^F(s) . (12) 

at asymptotically large values t. Let us find the limit distribution density 
foo(r) for the model dependence P(x), 

P{p/r) = c{^)\ a>0. (13) 

In this case, 

p(s) = — — , // = p(4) = 



a + s ' a + 4 

Thus, the average total surface area and the average size are defined by the 
formulas 

A(t) ^ (q + l)(a + 4) rl (15) 
and the average fragment number is 

Arm ^ * 3 . 

K) (a + 4) 3 (a+l)(a + 2)(a + 3) 
Using (13), let us transform the equation (12), 

F(s + 1) = ^rF(s) . 
a + 1 



4 



The solution of this functional equation which satisfies to the normalization 
condition F(l,t) = 1 is 

F(s) T(a + S) 



T{a+ l)(a+ I) 3 - 1 ' 

Computing the inverse Mellin transformation, we find the limit distribution 
on fragment size density 

(a + i)«+i r « e -(«+i)*- 
1 [a + 1J 

It is in accordance with the result of the work [3] which has been obtained in 
the case when a is the integer. However, the last distribution density differs 
essentially from the logarithmically normal Kolmogorov law. 
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